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We present Dalitz plot analyses for the decays of B mesons to D~D°K + 
and D°D°K + . [Charge conjugate reactions are implicitly assumed through¬ 
out.] We report the observation of the D* 1 (2700) + resonance in these two 
channels and obtain measurements of the mass M(D* 1 ( 2700) + ) = 2699hj. 
MeV/c 2 and of the width r(D* L (2700) + ) = 127f 2 g MeV, including statisti¬ 
cal and systematic uncertainties. In addition, we observe an enhancement 
in the D°K + invariant mass around 2350-2500 MeV/c 2 in both decays 
B° —» D~D°K + and B + —> D°D°K + , which we are not able to interpret. 
The results are based on 429 fb~ x of data containing 471 x 10 6 BB pairs 
collected at the T(4S) resonance with the BABAR detector at the SLAC 
National Accelerator Laboratory. 


PRESENTED AT 


The 7th International Workshop on Charm Physics 

(CHARM 2015) 

Detroit, MI, 18-22 May, 2015 


1 Work supported by the U.S. Department of Energy. 



1 Introduction 


We analyze B decays to DDK final states via Dalitz plots in order to measure the 
mass and width of the D*j( 2700) + . This is the first time that Dalitz plot analyses 
have been performed for B —>• DDK decays. 

The data analyzed here were recorded by the BABAR detector at the PEP-II 
asymmetric-energy e + e~ storage ring operating at the SLAC National Accelerator 
Laboratory. This analysis uses the complete BABAR data sample collected at the 
T(4S') resonance corresponding to an integrated luminosity of 429 fb _1 [lj. The 
BABAR detector is described in detail elsewhere [2]. Our Monte Carlo simulation 
uses EVTGEN [3] to model the kinematics of B mesons and JETSET [4] to model 
continuum processes, e + e _ —>■ qq (q = u, d, s, c ). The BABAR detector and its response 
to particle interactions are modeled using the GEANT4 [3] simulation package. 

2 Data set and selection 

The selection and reconstruction of B° —> D~ D°K + and B + —>■ D°D°K + , along with 
20 other B —> D^D^K modes, is described in Ref. pT]. We reconstruct D mesons 
in the modes D° —> A' _ 7t + , A" _ 7r + 7r°, A' - 7r + 7r - 7r + , and D + —>■ K~7t + 7t + . For the 
mode B + —> D D°K + , at least one of the D° mesons is required to decay to K~n + . 
We use the beam-energy-substituted mass ( tues ) an d A A, the difference between the 
reconstructed energy of the B candidate and the beam energy in the e + e~ center-of- 
mass frame to assist in identifying signal and in case of multiple candidates per event 
we use the latter to make a selection. Finally, we keep only events with | AA| < 10 — 14 
MeV depending on the D final state pa]. 

We fit the tues distributions, as described in detail in Ref. [6] to obtain signal 
yields; however, since there are 22 B —s- D^D^K modes that can cross-feed each 
other an iterative procedure is employed to obtain 635 ±47 and 901 ±54 signal events 
for B° —> D~ D a K + and B + —> D D°K + , respectively [6jj. For the Dalitz analyses we 
use the 5.275 < ??tes < 5.284 GeV/c 2 region to obtain a total of 1470 events with a 
signal purity of (38.6 ± 2.8 ± 2.1)% for B° —> D~D°K + aad 1894 events with a signal 
purity of (41.6 ± 2.5 ± 3.1)% for B + —> D D°K + , where the first uncertainties are 
statistical and the second systematic. 

3 Dalitz plot analyses 

We use an isobar model formalism to perform the Dalitz plot analysis [7]. The 
formalism is well-known and we merely note here that the dynamical amplitude, a 
Breit-Wigner form, contains a multiplicative factor F r which is the Blatt-Weisskopf 
damping factor for the resonance. [8] 
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Figure 1: Fits of the tues data distributions [6j for the modes B° —>■ D A>°A' + (left) 
and B+ ->• D°D°K+ (right). 



Figure 2: Dalitz plots for the modes B° —>■ D A>°A' + (left) and B + —> 
D°D°K + (right). 


We extract the complex amplitudes present in the data (from their moduli pi and 
phases 0j), and the mass and width of the D* 1 (2700) + resonance using an unbinned 
maximum likelihood fit where the minimized negative twice-log-likclihood is called 
T and the A* 1 (2700) + resonance is the reference amplitude. We compare fits using 


A T = T-T X 


nominal • 


The initial values of the parameters are randomized and we perform 250 such 
fits to ensure stable convergence to a global minimum. Both the efficiency and the 
background are parameterized by binned values and we perform a 2-dimensional 
interpolation to obtain values at desired locations on the plot. 

The Dalitz plots for B° — > D~D°K + and B + — > D°D°K + a,re shown in Fig. [ 2 ] The 
known amplitudes that could give a contribution in the Dalitz plot for both modes 
are nonresonant events, the A>* 1 (2700) + meson, and the A>* 2 (2573) + meson, which can 
decay to D°K + , but has not been observed in B —> D^DK decays. The D*j(2860) + 
state is not included in the nominal fit. For the mode B + —>- D D°K + , the additional 
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contributions from charmonium states included in the fit are the "0(3770) meson, and 
the 0(4160) meson. 

In the following fits, the masses and widths of these resonances are fixed to their 
world averages [9], except for the 0* 1 (2700) + where the parameters are free to vary. 
The spin of this resonance is assumed to be 1. 

Preliminary fits with the components mentioned above fail to give satisfactory x 2 l n doi- 
Also, the low mass region between 2350 and 2500 MeV/c 2 is not well described, es¬ 
pecially for B + —> D°D°K + . 



Figure 3: Projections of the Dalitz plot on the D°K + axis for the data (dots) and for 
the result of the preliminary fit (total histogram) for the modes B° —> D~D°K + ( left) 
and B + D°D°K+( right). 


After verifying that this enhancement is due to signal and not background or a 
cross-feed reflection, using two methods to subtract background, we proceed to fit 
it using an ad-hoc function that exponentially falls with m 2 (D°K + ). We call these 
our nominal fits. Since this region does not overlap significantly with the 0* 1 (2700) + 
region of interest, we are content to ascribe a systematic error to this effect. [We 
tried a low mass resonance but it only marginally improves the nominal fit.] The 
nominal fit for B° —> D~D°K + returns X 2 / n doi = 56/45 and the nominal fit for 
B + —>• O°O 0 A' + gives X 2 / n dof = 86/48. These fits are presented in Figures 4, 5 J and 
[6} The high value of the x 2 / n doi can be partly explained by differences in the data 
and fit densities in the 0(3770) region. 

We consider several sources of systematic uncertainties in the fit parameters such 
as the moduli, the phases, the fit fractions, and the mass and width of the i7* 1 (2700) + ; 
the combined result of this investigation is reported as the second error below. 


4 Results 


The results for the Dalitz plot analysis of the modes B° —>■ D O 0 A" + and B + —> 


D°D°K + ave presented in Tables [l] and 2 This is the first time the .0*027000 


is 
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Figure 4: Projections of the Dalitz plot on M(DD) axis for the data (dots) and for 
the result of the nominal fit (total histogram) for the modes B° —> D~D°K + ( left) 
and B + -»■ D°D°K+ (right). 



Figure 5: Projections of the Dalitz plot on M(DK) axis for the data (dots) and for 
the result of the nominal fit (total histogram) for the modes B° —> D~D°K + (left) 
and B + -»■ D°D°K+( right). 


observed in the decay B° —y D~D°K + . An excess at low D°K + invariant mass is 
evident but we have been unable to determine its origin. 

Using the Dalitz fit fractions and the total branching fractions measured in a pre¬ 
vious publication [6j with the exact same data sample we obtain the results presented 
in Table |3l 

We list results for the mass and width of the D* 1 (2700) + meson in Table |4j Com¬ 
bining modes we obtain 

M(D* 1 (2700) + ) = 2699ly 4 MeV/c 2 , (1) 

r(D* 1 (2700) + ) = 127+^MeV, 

compatible with the world averages. Repeating our fits with the J = 0 and J = 2 
hypotheses (see results in Table [5]) we conclude that J = 1 is strongly favored; further 
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Figure 6: Projections of the Dalitz plot on M(DK ) axis for the data (dots) and for 
the result of the nominal fit (total histogram) for the modes B° —y D~D°K + ( left) 
and B + —* D°D°K + (right). 

Table 1: Results from the Dalitz plot fit (moduli, phases, and fractions) for B° —>■ 
D~D°K + . The different contributions are listed: the D* 1 (2700) + and D* 2 (2573) + 
resonances, the nonresonant amplitude and the low-mass excess described by an ex¬ 
ponential. The first uncertainties are statistical and the second systematic. 


Contribution 

Modulus 

Phase (°) 

Fraction (%) 

.D* 1 (2700) + 

1.00 

0 

66.7 ± 7.8±1;| 

DU 2573) + 

0.031 ±0.008 ±0.002 

277 ± 17±| 

3.2 ± 1.6+°i 

Nonresonant 

1.33 ± 0.63+^ 

287 ±21+}° 

10.9 ±6.6+^ 

Exponential 

6.94 ± 1.83±g;|§ 

269 ± 33i^ 

9.9 ± 2.91^3 

Sum 



90.6 ± 10.7ljj;£ 


assuming parity conservation in the resonant decays we deduce that the D* 1 (2700) + 
has J p = 1~. 

Finally, our fits do not favor inclusion of the D*j(2860) + and D*j(3040) + reso¬ 
nances in the final states B° —>■ D~D°K + a,nd B + —* D°D°K + . 
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Table 2: Results from the Dalitz plot fit (moduli, phases, and fractions) for 
B+ —y D°D°K + . The different contributions are listed: the £>*02700) + , £>* 2 (2573)+, 
■0(3770), and 0(4160) resonances, and the low-mass excess described by an exponen¬ 
tial. The first uncertainties are statistical and the second systematic. 


Contribution Modulus Phase (°) Fraction (%) 


£>*02700)+ 



1.00 



0 

38.3 

± 

5 0+°' 8 
J - u -6.2 

DU 2573)+ 

0.021 ± 

0. 

01 n+0-009 
uiu -0.003 

267 

± 

3O0L 7 

0.6 

± 

i.nst 

0(3770) 

1.40 

± 

0 21+ 0 ' 20 
U ’ Zi -0.24 

284 

± 

90+26 

z ^-30 

9.0 

± 

3-USt 

0(4160) 

0.78 

± 

0-201°;“ 

188 

± 

13^17 

6.4 

± 

3.iy» 

Exponential 

16.15 

± 

2.26+^74 

308 

±8lg 

44.5 

± 

6.2 t'd 

Sum 







98.9 

± 

9 2+ 2 ' 5 
y - z -7.o 
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Table 3: Summary of partial branching fractions. The first uncertainties are statistical 
and the second systematic. The notation B° —y D~D* 1 ( 2700) + [ D°K + ] refers to 
B° -)■ D~D* sl ( 2700)+ followed by ^(2700)+ ->• D°K+. 


Mode 


B (10 

- 4 ) 



B° -»■ D~D* al (2700)+ 

[D°/7+] 

7.14 

± 

0.96 

± 

0.69 

B+ -)■ D° J D* 1 (2700) + 

[D°/7+] 

5.02 

± 

0.71 

± 

0.93 

B° -)■ D-D* 2 (2573) + 

[D°/7 + ] 

0.34 

± 

0.17 

± 

0.05 

D 0 i4* 2 (2573) + 

[D°/7+] 

0.08 

± 

0.14 

± 

0.05 

-)• ^(3770)77+ [D 

°D°] 

1.18 

± 

0.41 

± 

0.15 

5+ ->■ ^(4160)77+ [D 

°D°] 

0.84 

± 

0.41 

± 

0.33 


Table 4: Mass and width of the i4* 1 (2700) + meson obtained from the Dalitz plot 
analyses of the modes B° —y D~D°K + and B + —> D D°K + . The first uncertainties 
are statistical and the second systematic. 


Mode 

Mass (MeV/c 2 ) 

Width (MeV) 

B° ->• D~D°K+ 

2694 ± 8^3 3 

145 ± 24l 22 

B+ D°D°K+ 

2707 ± 8 ± 8 

113 ± 2lt 2 6 


Table 5: Value of A T and x 2 / n dof for the hypotheses J = 0,1, 2 for the two modes. 
The nominal fit is presented in bold characters. 


Mode 

J = 0 

J= 1 

J = 2 


AT 7 x 2 /n do f 

xV^dof 

AJ 7 X 2 / n dof 

B° ->• D~D°K+ 

131 131/45 

56/45 

108 125/45 

B+ -)■ D°D°K+ 

63 137/48 

86/48 

99 145/48 
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